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Photon breeding in relativistic jets involves multiplication of high-energy photons propagating from the jet to the 
external environment and back with the conversion into electron-positron pairs. The exponential growth of the energy 
density of these photons is a super-critical process powered by the bulk energy of the jet. The efficient deceleration 
of the jet outer layers creates a structured jet morphology with the fast spine and slow sheath. In initially fast and 
high-power jets even the spine can be decelerated efficiently leading to very high radiative efficiencies of conversion of 
the jet bulk energy into radiation. The decelerating, structured jets have angular distribution of radiation significantly 
broader than that predicted by a simple blob model with a constant Lorentz factor. This reconciles the discrepancy 
between the high Doppler factors determined by the fits to the spectra of TeV blazars and the low apparent velocities 
observed at VLBI scales as well as the low jet Lorentz factors required by the observed statistics and luminosity ratio 
of Fanaroff-Riley I radio galaxies and BL Lac objects. Photon breeding produces a population of high-energy leptons 
in agreement with the constraints on the electron injection function required by spectral fits of the TeV blazars. 
Relativistic pairs created outside the jet and emitting gamma-rays by inverse Compton process might explain the 
relatively high level of the TeV emission from the misaligned jet in the radio galaxies. The mechanism reproduces basic 
spectral features observed in blazars including the blazar sequence (shift of the spectral peaks towards lower energies 
with increasing luminosity). The mechanism is very robust and can operate in various environments characterised by 
the high photon density. 
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1. Introduction 

A general consensus exists on the radiative processes 
responsible for the blazar emission: synchrotron, syn- 
chrotron self-Compton (SSC) and external radiation 
Compton (ERC) mechanisms. However, neither the 
distance from the central black hole to the gamma- 
ray emitting region, nor the mechanisms of the en- 
ergy dissipation and electron acceleration in rela- 
tivistic jets are well understood. The internal shocks 
model [1] as well as magnetic reconnection models [2] 
have been discussed as possible scenarios for the en- 
ergy dissipation mechanism. 

The arguments based on the gamma-ray trans- 
parency and the observed short time-scale variability 
constrain the gamma-ray emitting region in blazars 
to distances of about ~ 10 17 cm [3]. This distance is 
often associated with the broad-emission line region 
(BLR). In internal shock model, there is no physi- 
cally sound way of explaining this distance. In ad- 
dition, the internal shocks are rather inefficient ra- 
diators, unless huge fluctuations of the Lorentz fac- 



tors are involved [4]. Very rapid TeV variability de- 
tected from PKS 2155-304 [5] and Mrk 501 [6] re- 
quires small distances and/or small emission region 
size and makes the internal shock models question- 
able. The reconnection models are not yet at the level 
to predict the distance or efficiency. 

Recently, a novel photon breeding mechanism for 
dissipation of the bulk energy of relativistic jets was 
suggested [7, 8]. The mechanism involves multiplica- 
tion of high-energy photons propagating from the jet 
to the external environment and back with the con- 
version into electron-positron pairs. The exponential 
growth of the energy density in these photons is a 
super-critical process [9] powered by the jet bulk en- 
ergy. The general scheme of the converter mecha- 
nism, where conversion of charged particles to neu- 
tral ones and back plays an important role in par- 
ticle acceleration, was first discussed independently 
in [10, 11]. The numerical studies of the operation 
of the photon breeding mechanism in ultrarelativis- 
tic shocks [11], showed the high radiative efficiency 
of the bulk energy dissipation. A similar problem for 
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relativistic jets in active galactic nuclei was investi- 
gated by us [7, 8] using a ballistic model of the jet 
and a detailed treatment of particle propagation and 
interactions. We showed that a large fraction of the 
jet bulk energy can be converted into gamma-rays, if 
the following conditions are met: jet Lorentz factor 
Tj > 5-10, the presence of the ambient soft photon 
field and a transversal or a chaotic magnetic field in 
the jet and the ambient medium. 

In this paper we discuss the observational pre- 
dictions of the photon breeding mechanism and com- 
pare them to the available data on blazars and radio 
galaxies. In Section 2 we discuss the possible sites, 
where photon breeding can operate. In Section 3 we 
describe the character of the jet radiation in the pho- 
ton breeding regime, particularly, the angular distri- 
butions of the emission components. We also show 
that it naturally resolves the Doppler factor crisis. 
In Section 4 we compare numerical simulations to 
the observed broad-band blazar spectra as well as to 
the electron distributions obtained from the spectral 
fits. 



2. Photon breeding operation sites 

The operation of the photon breeding mechanism re- 
quires relativistic motion with at least Tj > 4 [7] 
and a source of sufficiently dense transverse photon 
field (e.g. BLR). The photon breeding thus allows 
various emission sites. The minimum requirement is 
that the optical depth for the high-energy photons to 
pair production on soft photons is not much smaller 
than unity across the jet. If the background of soft 
photons is given by the multi-colour accretion disc of 
luminosity L^, this condition can be written as [7] a 



r 77 (i?) = 60 



45 



i?17@r 



<m > i, 



(i) 



where 6 m ax = fc7max/Tt c 2 w 10 -5 is the maximum 
disc temperature, R is the distance to the emission 
site, 9 = Rj/R is opening angle of the jet. 

At distances comparable to the size of the accre- 
tion disc, its direct radiation can serve as a target 
for high-energy photons and can trigger the photon 
avalanche, if the jet is already accelerated there. The 
photon-photon opacity through the jet is very large 



at these distances and the jet deceleration and radia- 
tive efficiency might be rather small, because only 
a very thin boundary layer can participate in pho- 
ton breeding. The presence of the X-rays from the 
accretion disc extending up to <~ 100 keV can be 
important and can trigger a pair cascade at small 
distances [12]. 

At larger distances the disc radiation becomes 
more beamed along the jet and the disc photons do 
not interact anymore with the high-energy radiation 
produced in the jet. The BLR radiation, however, 
supplies enough soft photons to satisfy the conditions 
for photon breeding. 

At the parsec scale, if the jet is still relativis- 
tic enough, the process can operate on the IR ra- 
diation from the dusty torus. Assuming a black 
body dust emission with the temperature Odust = 



^Tdust/™ c 2 « 10 7 L^4 5 R P c^ 2 , the optical depth i 
about [8] 



is 



'77 
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a We use standard notations Q 
dimcnsionlcss variables. 



10 x Qx in cgs units and for 



where rjd is the ratio of the dust to the disc luminosi- 
ties. Thus, even for ijd of the order of a few percent, 
the photon breeding is effective at the parsec scale 
in bright quasars. 

Even further out, if jet is still highly relativistic, 
photon breeding might operate on the stellar radia- 
tion field at kpc scale and at ^100 kpc on the cosmic 
microwave background radiation. 

Finally, we would like to note, that the emission 
from the blazar gamma-ray zone is very bright within 
the jet opening angle and can strongly affect the con- 
ditions for photon breeding at larger distances, where 
the radiation of the produced pairs can be responsi- 
ble for the jet emission observed in radio galaxies. 

3. Angular distribution, Doppler factor 
crisis and volume dissipation 

3.1. Emission pattern in photon 
breeding model 

Photon breeding in the relativistic shear layers 
causes efficient momentum exchange between rapidly 
moving outer layers of the jet and the external en- 
vironment [7, 8]. ERC scattering by the relativis- 
tic electron-positron pairs born within the jet, be- 
ing rather anisotropic in the jet comoving frame, not 
only produces gamma-ray emission, but also carries 
away a fraction of the jet momentum. This causes 
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Fig. 1. Distribution of the fluid Lorentz factor in the direc- 
tion of jet propagation (measured in units of radius of the jet 
cross-section Rj = 9R) for various starting jet Lorentz fac- 
tors: Ij = 14, 20, 30, 40. Solid curves show the Lorentz factor 
at the jet axis r = 0, and the dashed curves are for the jet 
boundary at 0.99-Rj. Parameters of the simulations: disc lumi- 
nosity L<j = 3 X 10 44 erg s _1 , jet power Lj = L<j, Poynting flux 
Lb = 0.2L(j) distance to the BLR R = 10 17 cm, jet opening 
angle 6 = 0.05. 

first deceleration of the jet's outer layers (see dashed 
curves in Fig. 1). As the cascade develops within the 
jet, even the spine can then be decelerated (see solid 
curves in Fig. 1) 

The main radiative mechanisms responsible for 
the jet emission are the same as usually assumed to 
describe blazar spectra: synchrotron, SSC, and ERC. 
In a differentially decelerating jet, the synchrotron 
radiation from the slower sheath can dominate the 
energy density of soft photons within the faster spine. 
Thus the energy loss may be dominated by the ex- 
ternal synchrotron Compton (ESC) mechanism (see 
also [13, 14]). 

The angular distribution of radiation from the 
differentially decelerating jet is very different from 
that predicted by simplistic models invoking rela- 
tivistically moving blobs. For an isotropically emit- 
ting blob, the angular pattern is given by 

m = f , (3) 
1 j 

where 6 = 1/Ij(l — (3jCosd) is the Doppler factor 
(see solid curves in Fig. 2), while the ERC beaming 
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Fig. 2. Angular distribution of luminosity in the three energy 
bands (IR-optical, X-rays and GeV gamma-rays) from the 
jet differentially decelerating under action of photon breeding 
(histograms, arbitrary normalization). Parameters: Tj = 20, 
L d = Lj = L B = 3 x 10 45 erg s _1 , R = 10 17 cm, 9 = 0.05. 
Theoretical dependencies expected for the steady relativistic 
jet emitting isotropically in the comoving frame (3) and the 
ERC radiation (4) are shown by solid curves (cases with T = 2 
and 20) and dashed curve, respectively. 

is given by [8, 15, 16] (dotted curve in Fig. 2) 

m = \L (4) 

Photon breeding gives a much wider beams. The his- 
tograms in Fig. 2 show the angular dependence of the 
luminosity in three typical energy bands (IR-optical, 
X-rays and gamma-rays). 

We see that the distribution of high-energy pho- 
tons at small angles follows quite closely equation 
(3). The low-energy (IR-optical and the X-rays) pho- 
tons are much less beamed. The spatial gradients of 
r is the main cause of that. The soft band is domi- 
nated by the synchrotron radiation and the emission 
at 6 = 0.5 exceeds the simple estimate (3) by 3 orders 
of magnitude. 

As a result of photon breeding, a significant pop- 
ulation of relativistic pairs is also born in the external 
environment. They inverse Compton scatter the disc 
and BLR radiation and synchrotron photons from 
the jet more or less isotropically producing an addi- 
tional component clearly visible at large angles. 
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3.2. Doppler factor crisis in TeV 
blazars 

The Doppler factors required by the homogeneous 
SSC models to describe the spectra of the blazars 
emitting at TeV energies (to avoid absorption by the 
infrared radiation), are large r$^20-100 [29-32]. On 
the other hand, the apparent velocities observed at 
the parsec scale in TeV blazars [17-19] are mildly rel- 
ativistic. In addition, the jet Lorentz factors, derived 
by matching the luminosity functions and statistics 
of BL Lac objects and their parent population of FR I 
radio galaxies, was estimated to be Tj ~3-6 [20-22]. 
Such a clear disagreement between the jet Lorentz 
factors determined by various methods has prompted 
to talk about the Doppler factor crisis [23]. 

One proposed solution to this crisis is that the 
jet Lorentz factor drops from 20-50 at the gamma- 
ray emitting, subparsec scale to just a few at the 
radio-emitting VLBI parsec scale [13]. For the view- 
ing angle within the beaming pattern of the initial 
jet, 8 < 1/Tj, the apparent velocity of the decelerated 
jet with Lorentz factor Tf <C Tj is [8] 

/3 f , app « 20r f 2 < 2^. (5) 
1 j 

Thus, deceleration of the jet to Tf < y/Tj/2 guar- 
anties that the apparent motion is subluminal. 

It was also suggested that the jet may consists of 
a fast spine and a slow sheath (so called structured 
jet) [14, 21]. In this model a slower sheath dominates 
the emission in the off-axis sources. 

However, both proposals lack a physical mech- 
anism for dissipation of a significant fraction of the 
jet bulk energy to achieve a required deceleration. 
The photon breeding mechanism, on the other hand, 
produces a decelerating and structured jet in a self- 
consistent way (Fig. 1; [8, 9]), and therefore it is ca- 
pable of unifying BL Lacs with radio galaxies keeping 
high Tj to explain the gamma-ray emission and may 
resolve thus the Doppler factor crisis. 

3.3. Off-axis emission and the TeV 
emission from radio galaxies 

As we have discussed above, the photon breeding 
mechanism naturally produces very broad beams of 
photons (see also [10, 24]). Interestingly, the ratio 
of soft to hard photon luminosities in not a mono- 
tonic function of the viewing angle (see Fig. 2). In 



slightly misaligned objects the gamma luminosity 
may be very weak compared to the optical luminos- 
ity, while again in radio galaxies observed at large 
angles ~ 1, the relative contribution of the gamma- 
ray flux increases due to the emission by the exter- 
nally produced pairs. The luminosity ratio between 
the nearly isotropic emission and the beamed emis- 
sion at 9 « 1/Tj is about Tr 4 (two powers of Tj 
come from the energy amplification in the jet and two 
powers appear because of beaming into a solid angle 
- 1/r?). This ratio is - T?/10 larger than that pre- 
dicted by equation (3). Thanks to the emission from 
the external medium, the nearby misaligned jets be- 
come observable in high-energy gamma-rays. 

The central source of M87, which has the best 
studied jet, has a rather low disc luminosity of only 
~ 10 42 ergs _1 [25]. Sufficient pair-production opacity 
(1) can be achieved at R < 10 16 cm, which is only ~ 
10 Schwarzschild radii for the estimated mass of the 
central black hole ~ 3T0 9 M Q [26]. At such a distance 
the photon breeding can operate on the direct disc 
radiation. The implied distance is consistent with the 
detection of the rapid (~ 10 5 s) variability of the TeV 
photon flux from M87 [27]. 

3.4. Observational appearance and 
volume dissipation 

Most of the emission from the jet undergoing photon 
breeding comes from the regions of maximum gradi- 
ent of r. In the steady-state, the volume emissivity 
contributing to the observed emission (at small an- 
gles to the jet) is 



In the case of small luminosities (and small Tj) only 
the outer layers of the jet suffer significant decelera- 
tion, resulting in the limb-brightening, with most of 
the emission coming from large z, where the cascade 
had time to develop fully. The situation is more com- 
plicated at high luminosities (or large Tj). The outer 
layers decelerate rapidly (see curves for Tj = 30, 40 
in Fig. 1) making the emission limb-brightened at 
small z, while the layers close to the jet axis deceler- 
ate later. Thus the emission from the core dominates 
at large z. At angles larger than > 1/rj, only slower, 
significantly decelerated layers contribute to the ob- 
served luminosity, and therefore we always would see 
the limb-brightened emission. 



16:26 



WSPC - Proceedings Trim Size: llin x 8.5in poutanen 



Our predictions are consistent with the observed 
limb-brightened morphology of the Mrk 501 jet at the 
parsec scale [28] . Such a structure could a result from 
the efficient deceleration and loading by relativistic 
pairs of the jet outer layers by the photon breeding 
at subparsec scale. 

In spite of the fact, that deceleration is more 
efficient at the jet boundary, the energy dissipation 
can involve the whole volume of the jet as it is defined 
by the mean-free path of high-energy photons. This 
is in a striking contrast to the Fermi-type diffusive 
acceleration, where only a narrow boundary layer is 
active, because of the much smaller Larmor radii. 

4. Spectral energy distribution of 

blazars and the electron distribution 

As a result of photon breeding, a population of high- 
energy electron-positron pairs is produced. The pair 
injection spectrum should be bounded and mirror 
(relative to m c c 2 ) the spectrum of the soft photon 
background, because the pairs in the jet are produced 
only by those high-energy photons that can interact 
with the soft photons. For example, if soft photons 
are from the accretion disc of maximum temperature 
©max, the low-energy injection cutoff in the jet frame 
is at 7 m ; n ~ Fj/60 max [7, 8]. This cutoff should de- 
pend weakly on the luminosity and the black hole 
mass, which defines the characteristic emission ra- 
dius, so that 6 max fa 10- 5 L^ 4 4 5 (M BH /10 8 M Q )- 1 /2 
and thus we can predict that 

- ~ 3 -° 5 (S)^ 4 fef ■ m 

At high luminosities, however, the synchrotron pho- 
tons from the jet can provide enough opacity and 
7min may be smaller. 

The detailed fitting of the broad-band spectra 
of TeV blazars and low-power BL Lacs with the SSC 
model shows that the distribution of relativistic elec- 
trons injected to the system can be modelled as a 
power-law of index p w 2-2.5 with the low-energy 
cutoff at 7 m ; n ~ 10 4 -10 5 and the high-energy cutoff 
at 7 max - 10 6 -10 7 [29-32]. b Such a peaked distribu- 
tion of injected electrons is difficult, if not impossible, 



b The high-luminosity objects require much lower electron en- 
ergies [29], but this result may be biased because the electron 
distribution is obtained from one-zone models, while the syn- 
chrotron peak in reality may not be related to the gamma-rays. 
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Fig. 3. Electron (and positron) energy distributions in the jet 
comoving frame (averaged over the jet volume) as a result of 
photon breeding. Simulation correspond to Tj = 20, Lj = L d , 
Lb = 0.2L(j. with L<j varying (from bottom to top histograms) 
from 5 X 10 43 to 10 46 erg s -1 . 

to obtain with the standard Fermi-type acceleration 
mechanisms, while the photon breeding can repro- 
duce it easily (see Fig. 3). 

At higher luminosities (and compactnesses), the 
high-energy photons produce a pair cascade. This 
reduces the mean energy per pair (see Fig. 3). At 
these compactnesses the pairs also cool more effi- 
ciently producing a cooling distribution dN/dj oc 
7~ 2 , which changes to a steeper j~ 3 behaviour at 
higher energies as a result of pair cascade [33] . The 
latter part of the pair energy distribution is responsi- 
ble for the flat part of photon spectrum (in E 2 dN/dE 
units). The pair cascade does not degrade the pair 
energies infinitely. When the energy of the photon in 
the jet frame drops below l/0 max Fj ~ 10 4 , it can 
easily escape and the cascade stops. The reduction 
of the mean pair energy results in an increased role 
of ERC, because lower-energy pairs interact with ex- 
ternal radiation in Thomson regime. 

The photon spectra predicted by the photon 
breeding model [7-9] demonstrate two-component 
structure (synchrotron plus Compton). The peaks 
shift to lower energies as luminosity increases (so 
called blazar sequence, [34]). At low luminosities the 
mean Lorentz factor of the electron-positron distri- 
bution is large (see Fig. 3), and SSC dominates over 
ERC because of the Klein-Nishina effect. For higher 
luminosities, ERC dominates the cooling. This be- 
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haviour is in general agreement with the observa- 
tional trends. However, there are some discrepancies 
between the shapes of the theoretical and the ob- 
served spectra. In our model at low luminosities, the 
synchrotron cutoff in the 100 MeV-1 GeV range is 
produced by very energetic pairs which are created 
in the jet as a result of interaction with very soft pho- 
tons. It is possible that the BLR spectrum is rather 
peaked, then the injection spectrum of pairs is nar- 
rower and softer, reducing thus the energy of the syn- 
chrotron cutoff. 

At high luminosities, our model predicts a cool- 
ing spectrum peaking around 100 MeV, while the 
blazars associated with quasars have two, clearly 
separated components. However, there is no con- 
sensus whether the two components are co-spatial. 
The radio-to-optical emission in the high-luminosity 
sources can be produced is a separate region, far 
from the gamma-ray emitting zone. It is also possi- 
ble that the pairs produced in the jet can be heated 
at larger distances producing the low-energy syn- 
chrotron component. 

5. Conclusions 

Photon breeding mechanism can operate at various 
sites characterised by sufficiently high photon den- 
sity such as in the vicinity of the accretion disc, the 
broad-emission line region, at parsec scale in the in- 
frared radiation field of hot dust, at kpc scale in the 
stellar radiation field and even at 100 kpc scale us- 
ing cosmic-microwave radiation as a target. As a re- 
sult of operation of photon breeding a decelerating, 
structured jets is produced. The angular distribution 
of radiation is much broader than that predicted by 
a blob model moving with a constant Lorentz factor. 
The broad emission pattern reconciles the discrep- 
ancy between the high Doppler factors determined 
by the fits to the spectra of TeV blazars and the low 
apparent velocities observed at VLBI scales as well 
as low jet Lorentz factors required by the statistics 
and luminosity ratio of FR I radio galaxies and BL 
Lac objects. The inverse Compton emission from the 
pairs produced outside the jet predicts a high level 
of gamma-ray emission from the misaligned jets in 
radio galaxies. 

Photon breeding effectively accelerates predom- 
inantly high-energy leptons in agreement with the 
constraints on the electron injection function deter- 
mined by the spectral fits to TeV blazars. The mecha- 



nism also reproduces basic spectral features observed 
in blazars including the blazar sequence. 
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